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ABSTRACT 

Aims. This paper presents a detailed and precise study of the characteristics of the Exoplanet Host Star and CoRoT main target HD 
52265, as derived from asteroseismic studies. The results are compared with previous estimates, with a comprehensive summary and 
discussion. 

Methods. The basic method is similar to that previously used by the Toulouse group for solar-type stars. Models are computed with 
various initial chemical compositions and the computed p-mode frequencies are compared with the observed ones. All models include 
atomic diffusion and the importance of radiative accelerations is discussed. Several tests are used, including the usual frequency 
combinations and the fits of the echelle diagrams. The possible surface effects are introduced and discussed. Automatic codes are also 
used to find the best model for this star (SEEK, AMP) and their results are compared with that obtained with the detailed method. 
Results. We find precise results for the mass, radius and age of this star, as well as its effective temperature and luminosity. We also 
give an estimate of the initial helium abundance. These results are important for the characterization of the star-planet system. 
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In the last few years the number of observed exoplanets have increased dramatically, owing to missions like CoRoT (Baglin et al. 
2006 ) and Kepler (Koch et al. 120101) . as well as many other ground based and space missions (see the exoplanet encyclopedia for 
a complete summary Q). Precise characterization of exoplanet host stars become more and more important in the framework of the 
detailed studies of the observed planetary systems. Constraints on the parameters and internal structure of the star can be obtained 
by comparing models with photometric and spectroscopic observations (Southworth l201 ll Basu et al. 120121 1. but the best precision 
is obtained from asteroseismology, when the stellar oscillations may be observed and analysed. This was the case, for example, for 
the exoplanet host stars (hereafter EHS) l Hor (Vauclair et al. 2008) and p Arae (Soriano & Vauclair [2010) . both observed with 
HARPS, as well as the EHS HAT-P-7, HAT-P-11 and TrES-2 (Christensen-Dalsgaard et al. l20T0b . observed with Kepler. 

Among the EHS, the star HD 52265 is one of the most precisely observed for asteroseismology, as it was the only EHS observed 
£3 | as a main target by CoRoT. This GOV metal-rich main sequence star has an orbiting jupiter-mass planet at 0.5 AU with a period of 
■ ■' 1 19 days (Naef et al. 1200 II Butler et al. l2000l l. It was continuously observed between December 13, 2008 and March 3, 2009, that 
is 1 17 consecutive days. As a result, 31 p-mode frequencies were reported, between 1500-2550 pHz, corresponding to t — 0, 1 and 
2 (Ballot et al. 1201 11 1. From this analysis, a large separation of < Av > = 98.4 + 0.1 pHz and a small separation of < 5vo2 >=8.1 
+ 0.2 pHz were found, and a complete asteroseismic analysis including mode lifetimes was presented. An extensive study of the 
seismic rotation of HD 52265 was performed by Stahn [20TTI and Gizon et al. 120 121 

Spectroscopic observations of this star were done by several groups, who gave different values of the observed triplet ([Fe/H], 
log g, Teff). Their results are given in Table 1. Some of these groups also observed lines of other elements, and gave detailed 
relative abundances. The results show an overall surmetallicity in this star, similar for most heavy elements, with a small dispersion. 
The Hipparcos parallax is 34.54 + 0.40 mas (van Leeuwen 2007), which leads to a luminosity value log L/Lq = 0.29 + 0.05. 
A spectroscopic follow up was also done with the Narval spectropolarimeter installed on the Bernard Lyot telescope at Pic du 
Midi Observatory (France) during December 2008 and January 2009, i.e., during CoRoT observations. No magnetic signature was 
observed. 



http://exoplanet.eu/index.php 
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Preliminary modeling of this star, using spectroscopic constraints, was done by Soriano et al. (2007 ), as a preparation to CoRoT 
observations. Evolutionary tracks were computed using the Toulouse-Geneva evolution code (TGEC). According to the spectro- 
scopic constraints, eight models with masses between 1.18 and 1.30 Mq and metallicities ranging from 0.19 to 0.27 were chosen 
for further analysis, and adiabatic p-modes frequencies were computed. Echelle diagrams for each selected model were presented 
as well as their corresponding large and small separations. A large separation around ~ 100 fiHz was predicted from these models 
except for the Takeda et al. (2005 ) values, which corresponded to a smaller large separation (around ~ 75 yuHz). 

The detailed CoRoT observations allow going further in this analysis, using the precise seismic results. First of all, the Takeda 
et al. (120051) values,which correspond to a more evolved star, are excluded. Now we present a complete asteroseismic analysis for 
HD 52265, and give precise results on the stellar parameters. 

The method and models used for the asteroseismic comparisons with observations are described in section 2. All the models, 
as described in section 2.1, include element gravitational settling. The seismic tests are discussed in section 2.2 and we discuss the 
influence of radiative accelerations on heavy elements in section 2.3. The results are given in section 3. Section 3. 1 is devoted to the 
results obtained without taking surface effects into account. We test the influence of varying the initial chemical composition and the 
introduction or not of radiative levitation on heavy elements. An analysis of surface effects and the consequences on the results is 
given in section 3.2. A first discussion of the results is given in section 3.3. Finally, in section 4 the results obtained using automatic 
codes to find the best models for this star from seismology (SEEK and AMP) are presented and compared with the previously 
obtained solutions. Summary and discussion are given in section 5. 

Table 1. Summary of previous spectroscopic studies of HD 52265 



[Fe/H] 


Trff 




l°g£ 


Reference 


0.27 


±0.02 


6162 ± 


22 


4.29 


±0.04 


Gonzalez et al. 2001 


0.23 


±0.05 


6103 ± 


52 


4.28 


± 0.12 


Santos et al. 2004 


0.19 


±0.03 


6069 ± 


15 


4.12 


±0.09 


Takeda et al. 2005 


0.19 


±0.03 


6076 ± 


44 


4.26 


±0.06 


Fischer & Valenti 2005 


0.24 


±0.02 


6179 ± 


18 


4.36 


±0.03 


Gillon & Magain 2006 


0.19 


±0.05 


6100 ± 


60 


4.35 


±0.09 


Ballot etal. 2011 



2. Computations with TGEC 

2.1. Stellar models 

Stellar models were computed using the TGEC code, (Hui-Bon-Hoa l2008l Theado et al. I2012I ). with the OPAL equation of state 
and opacities (Rogers & Navfonov I2002I Iglesias & Rogers 1996), and the NACRE nuclear reaction rates (Angulo et al. [1999 1. 
Convection was treated using the mixing length theory. For all models, the mixing length parameter was adjusted to that of the solar 
case, i.e. a = 1.8 without overshooting and without extra-mixing. Gravitational settling of helium and metals was included using 
the Paquette prescription (Paquette et al. [1986; Michaud et al. 2004). Radiative accelerations of metals were also introduced using 
the SVP method (Single Valued Parameters approximation, see Alecian & LeBlanc 120021 LeBlanc & Alecian 2004 and Theado et 
al. 120091 ). As most stellar evolution codes neglect these radiative accelerations, we analysed the effects on the seismic results of 
introducing them or not. 

Evolutionary tracks were computed for two metallicity values and two different initial helium abundances. The metallicity values 
were chosen as [Fe/H] = 0.23 and 0.30, so that after diffusion the final model value lies inside the observed range. Here [Fe/H] 
represents the global overmetallicity with respect to the Sun, defined as [log(Z/X) t - log(Z/X)o], where Z and X are computed 
at the stellar surface. Considering the small dispersion of the detailed abundances, this value may be compared with the observed 
[Fe/H]. A discussion of the computed detailed abundance variations is given in section 2.3. The initial helium values are labeled as 
Yq and Yg, where Y0 is the solar helium value taken from Grevesse & Noels [T9"9"3"l and Yg is a helium abundance which increases 
with Z as expected if the local medium follows the general trend observed for the chemical evolution of galaxies (C.f. Izotov & 
Thuan l2004l and l20T0l ). 

Adiabatic oscillation frequencies were computed for many models along the evolutionary tracks using the PULSE code 
(Brassard & Charpinet 2008 I. We computed these frequencies for degrees 6 = to i = 3. For comparisons with the observations 
(seismic tests) we used the same frequency interval for the computed as for the observed frequencies for consistency, i.e. between 
1500 and 2550 yuHz, as discussed below. 

2.2. Seismic Tests 

A well known characteristic of p modes in spherical stars is that modes of the same degree with successive radial order n are nearly 
equally spaced in frequency (e.g. Tassoul 1980). The large separation is defined as: 

AV„^ = V n+ l/ - V n / (1) 

In real stars, this large separation slightly varies with frequency, so that an average value has to be used for comparisons between 
models and observations. One has to be careful to use the same frequency range in both cases to do the comparisons. Taking this 
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Fig. 1. Evolutionary tracks in the log g vs log plane for the various sets of metallicity and helium abundance, for a — 1.8 (see 
text for details). The symbols indicate the error boxes of Gonzalez et al. (2001 ) (asterisks), Santos et al. (2004 ) (diamonds), Gillon 
& Magain (2006) (squares), Fisher & Valenti (2005 ) (triangles), and Ballot et al. (1201 II ) (crosses). The straight thick line represents 
the iso-< Av > line, with < Av >=98.4 //Hz. 



into account, the fit between the computed and measured large separations is the first step in the process of comparisons. The large 
separation gives access to the stellar average density (Urrich [T9861 White et al. 1201 II ). 

A second characteristic of p modes is that the difference between (n, () modes and (n — \,( + 2) ones varies very slowly with 
frequency. The small separations are defined as: 

Sv n<( = v n ,« - V n -i/+2 (2) 

These small separations are most sensitive to the stellar core and may give information on the extension of the convective core 
in some stars (Tassoul 1980; Roxburg & Vorontsov |1994| Gough 1986; Soriano & Vauclair2008). 

Provided that the stellar chemical composition is precisely known, the knowledge of both the large and the small separations, 
which may be plotted in the so-called C-D diagrams, gives good constraints on the stellar parameters (Christensen-Dalsgaard 1984, 
White et al. 1201 II ). However, whereas the stellar metallicity can be precisely derived from spectroscopy, the helium content of solar- 
type stars is not directly known from observations. This ignorance leads to important uncertainties on the evolutionary tracks, and 
thus on the derived stellar parameters. 

We analysed these uncertainties in detail by computing models with various chemical compositions. For each stellar evolutionary 
track that we computed, we first searched for the model which had an average large separation of < Av > = 98.4 + 0.1 //Hz. As the 
large separation continously decreases when the star evolves along the main sequence, one model only is found with the observed 
value (within the uncertainties). Then, for each set of computations done with a given initial chemical composition ([Fe/H] and Y), 
we looked for the model which best fitted the small separations observed between modes of I -2 and 0, < 6vq2 >=8.1 + 0.2 //Hz. 
We further proceed with detailed comparisons of observed and computed echelle diagrams. 

The final comparison between the models and the seismic observations needs taking into account the surface effects induced by 
the waves behavior in the outer stellar layers. We computed the frequency shift induced by such effects, using the recipe proposed 
by Kjeldsen et al. (2008). In this case, the large separations are modified, as discussed in section 3.2, which leads to corrections in 
the results. 
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Fig. 2. HR diagrams for a — 1 .8 (see text for details). Error boxes are the same horizontally as presented in FigureQ] for the effective 
temperatures, and correspond vertically to the luminosity uncertainty. The straight thick line represents the iso-< Av > line, with 
< Av >=98.4 fMz. 



2.3. Test on atomic diffusion: radiative accelerations 

Atomic diffusion is a very important process inside stars: it can modify the atmospheric abundances and also have strong implica- 
tions for the internal structure of stars (Richard, Michaud & Richer 2001 1 Michaud et al. 120041 Theado et al. 2009 ). At the present 
time, most stellar evolution codes include the computation of gravitational settling, but not the computation of radiative accelerations 
of heavy elements. These accelerations, which oppose gravitation, are negligible for the Sun but they rapidly become important for 
more massive stars (Michaud et al. |19761 l. The variations with depth of the radiative accelerations on specific elements can lead to 
their accumulation or depletion in various layers inside stars. The presence of a heavy iron layer above layers with smaller molecular 
weights creates an inverse /z-gradient, unstable towards thermohaline convection (Vauclair 2004). This induced mixing has also to 
be taken into account in computations of stellar modeling (Theado et al. 120091 ). An improved TGEC version including radiative 
accelerations on C, N, O, Ca and Fe was recently developped (Theado et al. 2012 ). This new version was used to compare the os- 
cillation frequencies computed with and without introducing the radiative accelerations in the models. We found that, for solar-type 
stars (masses less than 1.30 M0), the difference in the computed frequencies is small. Two models with the same mass of 1.28 
Mq and the same average large separations of 98.26 yuHz, one computed with the radiative accelerations and one neglecting them, 
present differences in the average small separations of the order of 0.01 //Hz. These differences between the two models decrease 
with decreasing stellar mass. We conclude that the radiative accelerations may be neglected in the following computations. 

The result that radiative accelerations have no important consequences in the present models is consistent with the fact that the 
detailed abundances do not show large relative variations for the observed elements (Section 1). Indeed, in case of gravitational 
settling, the abundances of heavy elements all decrease in a similar way, in spite of their different masses, due to the slowing down 
effect of the coulomb forces. The diffusion velocities vary typically as A/Z 2 where A is the mass number and Z the charge, which is 
similar for various elements in stellar conditions. This behavior would be modified if radiative accelerations were important. 
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Table 2. Examples of models with a—l.S, without surface effects. 



[Fe/H],- 


Yi 


M/Mq 


Age 


[Fe/H] s 


Y 5 


l°g<? 


log T eff 


log (L/L Q ) 


R/R© 


M/R J 


< Av > 


< <5V 02 > 








[Gyr] 






[K] 








[solar units] 


[/(Hz] 


IjuHz] 


n 93 

U.ZJ 


n 9Q3 

U.Z7J 


1 1 8 
1 . lo 


^ 689 
J.OOZ 


n i 6 
u. 10 


PI 9zLA 


A 967 


3 778 
J. 1 1 o 


n 3 1 n 

U. J 1U 


1 398 
1 . JZo 




Q8 ^ 1 


7 08 


93 

U.Z J 


9Q3 


1 90 

1 .Z.U 


3 904 


16 


n 946 


4 771 


3 789 

J. 1 oz 


n W) 


1 333 
1 .jjj 


s 1 

U.J 1 


08 ^8 


7 S4 

/ ,J4 


93 

U.Z J 


9Q3 

U.Z;7 J 


1 99 
1 ,zz 




U.1U 


948 


A 11% 

H-.Z 1 J 


3 787 
J. 1 o / 


3ss 


1 3zL1 
1 . jt 1 


s 1 

U.J 1 


08 "36 


7 8S 


93 

U.ZJ 


9Q3 

U.Z;7 J 


1 94 


9 416 


f) 1 S 


9zL9 


A 976 
h-.z / o 


3 701 


37S 

U. J 1 J 


1 3zL7 
1 . jt / 


S 1 

U.J 1 


08 "34 


8 


n 93 

U.ZJ 


n 971 

U.Z / 1 


1 99 
1 .zz 


^ 7S6 

J. / JO 


O 1 6 
U. 10 


H 998 

U.ZZo 


A 979 
M-.Z / Z 


3 776 


PI 3 1 9 
U. J 1Z 


1 3zL3 
1 . JH-J 


U.j 1 


08 33 

70.JJ 


7 1 1 


O 93 
U.ZJ 


971 
U.Z / 1 


1 9A 


^ 98^ 


O 1 6 
U. 10 


PI 998 
U.ZZo 


A 97S 


3 780 
j. / oU 


n 33A 

U. J JH- 


1 3AQ 


s 1 

U.J 1 


Q8 3 1 

70.J1 


7 sn 


93 


971 

U.Z / 1 


1 96 

1 .Z-U 


? 86^ 


n 16 

U. IO 


930 
u.zju 


4 978 


3 78S 


3ss 

U.JJJ 


1 3SS 

1 . j j j 


S 1 

U.J 1 


08 34 


7 08 

/ .70 


93 

U.ZJ 


971 

U.Z / 1 


1 98 


9 461 


16 

u. io 


997 
u.zz. / 


4 981 


3 78Q 


37S 

U. J 1 j 


1 361 


S 1 

U.J 1 


08 3^ 
70.JJ 


8 


n 3n 
U.jU 


U.jUj 


1 1(\ 
1 .ZU 


^ 90Q 


o 


n 9<7 


/I 97 1 
4.Z / 1 


3 77Q 


n 3 1 a 


1 333 


U.J 1 


QQ A 9 


7 


0.30 


0.303 


1.22 


2.820 


0.23 


0.258 


4.273 


3.783 


0.339 


1.341 


0.51 


98.35 


7.75 


0.30 


0.303 


1.24 


2.431 


0.23 


0.261 


4.276 


3.787 


0.361 


1.347 


0.51 


98.38 


8.11 


0.30 


0.303 


1.26 


2.072 


0.22 


0.254 


4.282 


3.792 


0.382 


1.354 


0.51 


98.38 


8.59 


0.30 


0.271 


1.24 


3.771 


0.23 


0.231 


4.274 


3.771 


0.299 


1.349 


0.51 


98.32 


7.11 


0.30 


0.271 


1.26 


3.293 


0.23 


0.231 


4.277 


3.776 


0.320 


1.356 


0.51 


98.37 


7.64 


0.30 


0.271 


1.28 


2.865 


0.23 


0.231 


4.280 


3.780 


0.341 


1.363 


0.51 


98.34 


7.87 


0.30 


0.271 


1.30 


2.476 


0.24 


0.233 


4.282 


3.784 


0.362 


1.369 


0.51 


98.39 


8.25 



Table 3. Examples of models with ff=1.8 , including surface effects. Here < Av* > represents the Kjeldsen et al. ([2008 (-corrected 
large separations. 



[Fe/H], 


Y, 


M/Mq 


Age 


[Fe/H] s 


Y s 


l°gg 


log r cff 


log (L/L e ) 


R/R© 


M/R j 


< Av* > 


< Av > 


< 6V02 > 








[Gyr] 






[K] 








[solar units] 


[j"Hz] 


[juHz] 


buHz] 


0.23 


0.293 


1.18 


3.548 


0.16 


0.247 


4.275 


3.779 


0.306 


1.315 


0.52 


99.67 


98.27 


7.27 


0.23 


0.293 


1.20 


3.069 


0.16 


0.247 


4.279 


3.784 


0.328 


1.321 


0.52 


99.68 


98.27 


7.83 


0.23 


0.293 


1.22 


2.670 


0.16 


0.250 


4.282 


3.788 


0.350 


1.327 


0.52 


99.86 


98.19 


8.11 


0.23 


0.293 


1.24 


2.251 


0.16 


0.244 


4.287 


3.792 


0.369 


1.331 


0.52 


100.08 


98.19 


8.47 


0.23 


0.271 


1.20 


4.130 


0.16 


0.228 


4.276 


3.773 


0.287 


1.325 


0.52 


99.51 


98.19 


6.85 


0.23 


0.271 


1.22 


3.622 


0.16 


0.228 


4.279 


3.777 


0.309 


1.331 


0.52 


99.56 


98.18 


7.26 


0.23 


0.271 


1.24 


3.149 


0.16 


0.229 


4.283 


3.781 


0.330 


1.337 


0.52 


99.61 


98.18 


7.69 


0.23 


0.271 


1.26 


2.730 


0.16 


0.232 


4.286 


3.786 


0.351 


1.343 


0.52 


99.62 


98.07 


8.11 


0.30 


0.303 


1.20 


3.089 


0.23 


0.258 


4.278 


3.780 


0.313 


1.322 


0.52 


99.64 


98.32 


7.70 


0.30 


0.303 


1.22 


2.685 


0.23 


0.260 


4.282 


3.784 


0.335 


1.328 


0.52 


99.72 


98.13 


7.94 


0.30 


0.303 


1.24 


2.296 


0.24 


0.263 


4.285 


3.788 


0.357 


1.334 


0.52 


99.78 


98.15 


8.34 


0.30 


0.303 


1.26 


1.907 


0.23 


0.257 


4.290 


3.793 


0.376 


1.337 


0.52 


100.23 


98.20 


8.84 


0.30 


0.271 


1.24 


3.637 


0.23 


0.231 


4.282 


3.773 


0.295 


1.338 


0.52 


99.56 


98.27 


7.30 


0.30 


0.271 


1.26 


3.173 


0.23 


0.232 


4.284 


3.777 


0.317 


1.345 


0.52 


99.50 


98.16 


7.64 


0.30 


0.271 


1.28 


2.730 


0.23 


0.233 


4.288 


3.781 


0.337 


1.351 


0.52 


99.64 


98.09 


8.08 


0.30 


0.271 


1.30 


2.326 


0.23 


0.236 


4.291 


3.785 


0.357 


1.355 


0.52 


99.80 


98.20 


8.52 



3. Results 

3.1. Computations without surface effects 

The computations of evolutionary tracks, with two initial metallicity values, [Fe/H] = 0.23 and 0.30, and two different helium 
abundances, Y0 equal to 0.271 and Yc respectively equal to 0.293 and 0.303 for the two metallicity values, lead to four different 
sets of tracks, each of them covering a mass range from 1.10 to 1.30 Mq. A few of these tracks are presented in Figures 1 and 2. 
Error boxes of five of the spectroscopic studies given in Table 1 are also drawn in these figures. 

As explained in previous sections, we found, along each evolutionary track, a model which has an average large separation 
consistent with the observed one, computed in the same frequency range, of ~ [1.5, 2.5] mHz. The location of all these models in 
the log g - log T e ff plane as well as in the HR diagrams are indicated in Figures [TJ and [2] with iso-< Av > 98.4 fiHz lines. For each 
case, we also computed the average small separation, 6v„fl2 (Table 2). We can see that for models with the same large separation, the 
small separation increases for increasing mass, so that in each case (i.e. for each set of chemical composition), there is a model that 
is consistent with both the large and small separations. However, when comparing the absolute model frequencies with the observed 
ones, we find that we must shift the computed frequencies by about 20 /iHz to obtain the best fit with the observed ones. This offset 
is attributed to surface effects. The echelle diagrams corresponding to these best models are given in Figure 3. 

3.2. Computations including surface effects 

We know that stellar modeling fails to represent properly the near-surface layers of the stars. As a consequence, there is a systematic 
offset between observed and computed frequencies. This offset is independent of the angular degree I and increases with frequency. 
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Fig. 3. Echelle diagrams for the best model found for each set of evolutionary tracks calculated with a — 1.8. In this set of models, 
the surface effects are not included. The model frequencies (crosses) are compared to the observed frequencies, represented as 
diamonds (^=0), asterisks (t=l), and triangles ((-2). To obtain these best fits we had to shift the model frequencies by respectively 
(top left, top right, bottom left, bottom right) 27, 25, 24, 22 fiHz (see text) 

This has been studied in the case of the Sun, and a similar offset is expected to occur in other stars. Using the Sun as a reference, 
Kjeldsen et al. (2008) suggested that for other stars, the near-surface correction on the frequencies may be approximated by: 



Vobs(n) - Vbest(n) = a 



v bs(n) 



(3) 



where v b s (n) are the observed i-0 frequencies with radial order n, Vbest(n) are the calculated frequencies for the best model which 
is the model that best describes the star but which still fails to model correctly the near-surface layers, and vo is a constant reference 
frequency, chosen to be that of the frequency at the maximum amplitude in the power spectrum. The parameter a may be derived as 
a function of the parameter b, which has to be adjusted to the solar case. We used this method to find the frequency corrections that 
we have to apply to our models, which lead to a new average large separation < Av* >, slightly larger than the observed one. We 
then proceed as before to derive models for the same sets of values of the helium abundance and metallicity, using a = 1 .8 (Table 3). 
Figure |4]presents echelle diagrams obtained with the new best models. In these graphs, the frequency corrections have been applied 
to the new model frequencies in order to compare them directly with the observations. 



3.3. Best models and discussion 

We give in Table [4] the parameters of the best models that we obtained for the four different sets of chemical composition. These 
models have been computed with a mixing length parameter of 1.8, the frequencies have been corrected for surface effects and 
the computed large and small separations are the closest to the observed ones in the sample. We also give the x 2 values for the 
comparisons of the three I — 0, 1,2 lines in the echelle diagrams. 

Figure [5] displays the large and small separations as a function of the frequency for the observations and the four best models. 
The pattern observed in the large separations are well reproduced by the models. For the small separations, the agreement is also 
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Fig. 4. Echelle diagrams for the best model, including near-surface corrections, as proposed by Kjeldsen et al. 2008 found for each 
set of evolutionary tracks calculated with a — 1.8, in comparison with observed frequencies. The symbols are the same ones as in 
Figure 3. 



Table 4. Best models obtained with the TGEC code, including near surface corrections. 



M/Mq 


L/Lq 


R/R 


l°g£ 


Teff [K] 


age [Gyr] 


[Fe/HL 


Y, 


[Fe/H]. v 


Y s 


x 1 


1.22 


2.239 


1.327 


4.282 


6143 


2.670 


0.23 


0.293 


0.16 


0.250 


5.06 


1.23 


2.219 


1.330 


4.283 


6120 


2.476 


0.30 


0.303 


0.23 


0.262 


3.34 


1.26 


2.244 


1.343 


4.286 


6109 


2.730 


0.23 


0.271 


0.16 


0.232 


6.42 


1.28 


2.173 


1.351 


4.288 


6043 


2.730 


0.30 


0.271 


0.23 


0.233 


3.51 



very good except for two points at large frequencies (2284 yuHz and 2479 yuHz). This suggests that the uncertainties given in Ballot 
et al. (1201 11 1 for these points were underestimated. 

Several concluding points can already be derived from a first analysis of Table 4. First of all, the stellar gravity is obtained, as 
usual, with a precision of order 0.1%. The mass and age depend basically on the chosen value for the initial helium content. For a 
low helium value, the mass is 1.26 to 1.28 Mq, and the age 2.73 Gyr, whereas for a larger helium value the mass is slightly smaller 
(around 1.22 - 1.23 Mq) as well as the age (2.48 to 2.68 Gyr). In any case the radius and luminosity are known with a precision of 
order 1%. 

We can go further by comparing the effective temperatures of the models with the spectroscopic observations (Figure 6). As 
found before for the cases of l Hor (Vauclair et al. |?008 1 and fi Arae (Soriano & Vauclair 2010 1, the effective temperatures of the 
best models are smaller for smaller initial helium and larger for smaller metallicity. In the present case, we find that the model with 
the largest metallicity and the smallest helium, represented by black square in figure 6, is at the coolest limit of the observational 
boxes and thus may be excluded from the sample on spectroscopic considerations. We derive the stellar parameters from a mean 
value of the resulting three models, as given in Table 5. 
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Fig. 5. Comparisons between the large separations (top) and the small separations (bottom) of the four best models indicated by 
squares for models with [Fe/H] s =0.23 ([Fe/H],=0.30) and triangles for models with [Fe/H] s =0.16 ([Fe/H],=0.23). Empty symbols 
are used to indicated an helium abundance according to Y^and filled symbols for a solar helium abundance (see text for details.) 
Observations are represented as white diamonds. 



4. Scaling relations and automatic fits 

4.1. Scaling relations 

The empirical scaling relations proposed by Kjeldsen & Bedding (1995) can give approximate values of the mass and radius of a 
star from the observed average large separation, the frequency at the maximum of the power spectrum and the observed effective 
temperature: 



M _ / 135////z \ 4 / 
Mq \< Av > j \: 



3050fiHz) \5111K 



eff 



3/2 



R_ = 1 135 fiHz f I v m „ x \ I Teff \ 
R Q \< Ay > J \3050fiHzJ\5mK) 



1/2 



(4) 



(5) 



For HD 52265, the frequency at the maximum amplitude is v nmx = 2090 + 20 fjHz (Ballot et al. 201 1). With a large separation 
of 98.4 //Hz and an effective temperature of 6100 K, we obtain from these relations M = 1.23Mq and R = 1.32/?q. With a large 
separation corrected for the surface effects, of 99.4 /iHz, we obtain M - 1.19M0 and R = 1.29/?q. In spite of the uncertainties, 
these results are in good agreement with our own results. 



4.2. Results from the SEEK code 

Computations have been done for this star using the SEEK automatic code (Quirion et al. 120101 and Gizon et al. 12012) . This code 
makes use of a large grid of stellar models computed with the Aarhus Stellar Evolution Code (ASTEC). It searches for the best 
model corresponding to a seismically observed star, with the help of Bayesian statistics. The input parameters are the large and small 
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Fig. 6. Location of the best models in the log g vs log T e g plane. Triangles indicate models with [Fe/H] s =0.16 ([Fe/H],=0.23), and 
squares indicate models with [Fe/H] s =0.23 ([Fe/H],=0.30). Filled symbols are used to show the models with a solar helium value. 
Error boxes correspond to Gonzalez et al. (2001) (asterisks), Santos et al. ((2004) (diamonds), Gillon & Magain (2006) (squares), 
Fisher & Valenti (2005) (triangles), and Ballot et al. (1201 11 1 (crosses) spectroscopic studies. 



average separations, the spectroscopic observables (Jeff, log g, [Fe/H]) and the absolute magnitude. The output gives the stellar 
mass, the radius and the age. In the case of HD 52265 the values given by the SEEK code for the mass and radius are slightly larger 
than our results, and the age is smaller : M = 1.27 ± 0.03M Q , R = 1.34 ± Q.Q2Rq and age = 2.37 ± 0.29 Gyr. The difference between 
the SEEK results and ours may be related to a different initial helium, to slightly different values of the average large and small 
separations as given by Gizon et al. (120121 1. or to the fact that the SEEK results correspond to a secondary maximum of probability, 
as discussed below. 



4.3. Results from the Asteroseismic Modeling Portal 

We also performed computations for HD 52265 using the Asteroseismic Modeling Portal (https://amp.ucar.edu/). The AMP provides 
a web-based interface for deriving stellar parameters for Sun-like stars from asteroseismic data. AMP was developed at the High 
Altitude Observatory and the Computational & Information Systems Laboratory of the National Center for Atmospheric Research 
(Metcalfe et al. 120091 ) . It uses the ASTEC and ADIPLS codes (Christensen-Dalsgaard 2008 a,b) coupled with a parallel genetic 
algorithm (Metcalfe & Charbonneau 2003i. Two different computations were done, the first one, AMP(a), by S. Vauclair using all 
the observed frequencies as given by Ballot et al. (1201 1) , and the second one, AMP(b), by S. Mathur using only the most reliable 
frequencies ((1=0, n=14), (1=1, n=14) and (l=2,n=15) were excluded). The final results are very close to the parameters found by 
using the TGEC code (Table 6). Interestingly enough, the code found also solutions for a mass of 1.27 Mq with a small Y (about 
0.26) and a small age (about 2.7 Gyr) but the^ 2 tests showed that these results corresponded to secondary maxima, not to the best 
solution. 



Table 5. Final results for the parameters of the exoplanet-host star HD 52265 obtained with the TGEC code. 



M/M e = 1.24 + 0.02 
R/R = 1.33 + 0.02 
L/L = 2.23 ± 0.03 
log g = 4.284 ± 0.002 
Age (Gyr) = 2.6 + 0.2 



[Fe/Hl = 0.27 ± 0.04 
Y,- = 0.28 ± 0.02 

[Fe/H] s = 0.20 ± 0.04 
Y, = 0.25 ± 0.02 

T eff (K) = 6120 + 20 



5. Summary and Conclusions 

In the present paper, we performed a detailed analysis of the exoplanet-host star HD 52265, which has been observed by CoRoT 
during 117 consecutive days, as one of the main targets. The beautiful observational results obtained for this star (Ballot et al. 
1201 11 1 allowed a precise determination of its parameters, using classical comparisons between models computed with the TGEC 
and observational data. In our computations, we included atomic diffusion of helium and heavy elements. We found that, for the 
computed stellar models, the effects of radiative accelerations on individual elements is small and may be neglected. This result is 
consistent with the fact that detailed abundance analysis show similar enhancements in heavy elements compared to the Sun. We 
iterated the model computations so as to find a final surface [Fe/H],s in the observational range. We also compared these results 
with those obtained using approximate scaling relations (Kjeldsen & Bedding 1995), and automatic codes like SEEK (Quirion et al. 
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Table 6. Final results for the parameters of the exoplanet-host star HD 52265 obtained from AMP automatic analysis, using all 
observed frequencies (a) or only the most reliable frequencies (b). 



2010) and AMP (Metcalfe et al. 2009). Although the detailed physics included in the models is different, these results are in good 
agreement. 

The good concordance between the results obtained with the TGEC code and the AMP for Sun-like stars was already proved 
with the star fi Arae. The results for this star, which were published separately (Soriano & Vauclair 2010 for TGEC; Dogan et al. 
120 121 for AMP), are presented in the appendix for comparison. Altogether, these works represent real success of the asteroseismic 
studies for Sun-like stars. 

Acknowledgements. Much thanks are due to Travis Metcalfe for fruitful discussions and for introducing SV to AMP. Computational resources were provided by 
TeraGrid allocation TG-AST090107 through the Asteroseismic Modeling Portal (AMP). NCAR is partially supported by the Natiional Science Foundation. We also 
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Appendix A: 

Comparison of the results given by the TGEC analysis (Soriano & Vauclair 2008 ) and the Asteroseismic Modeling Portal (Dogan 
et al l2012l l for the exoplanet-host sun-like star fj. Arae. 



Table A.l. Comparison of the results obtained with the TGEC and the Asteroseismic Modeling Portail for fi Arae. 



TGEC 



AMP 



M/Mq 
R/Rq 
L/L 

logg 



1.1 +0.02 
1.36 ±0.06 
1.90 ±0.10 

4.215 ±0.005 
5820 ± 50 

+0.32 ± 0.02 
0.30 ± 0.01 
6.34 ± 0.80 



1.1 
1.365 
1.894 
4.213 
5807 
+0.418 
0.316 
6.66 



[Fe/H] 




Y 



age(Gyr) 



11 



